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Molecular spin crossover switches are the objects of intense theoretical and experimental studies 
in recent years. This interest is due to the fact that these systems allow one to control their spin 
state by applying an external photo-, thermo-, piezo-, or magnetic stimuli. The greatest amount of 
research is currently devoted to the study the effect of the photoexcitation on the bi-stable states 
of spin crossover single molecular magnets (SMMs). The main limitation of photo-induced bi-stable 
states is their short lifetime. In this paper we present the results of a study of the spin dynamics of 
the CoOEP molecule in the Low Spin (LS) state and the High Spin (HS) state induced by applying 
the magnetic pulse of 36.ST. We show that the spin switching in case of the HS state of the CoOEP 
molecule is characterized by a long lifetime and is dependent on the magnitude and duration of the 
applied field. Thus, after applying an external stimuli the system in the LS state after the spin 
switching reverts to its ground state, whereas the system in the HS state remains in the excited 
state for a long time. We found that the temperature dependency of magnetic susceptibility shows 
an abrupt thermal spin transition between two spin states at 40K. The proposed here theoretical 
approach opens the way to create modern devices for spintronics with the controllable spin switching 
process. 

PACS numbers: 


I. INTRODUCTION 

Spin engineering mostly concentrates on the manipu¬ 
lation of spin degrees of freedom in various magnetic ma¬ 
terials and the exploration of their potential applications 
in spintronics^. Current examples for the applications of 
spin engineering in magnetic nanostructures and devices 
are increasingly ample, ranging from the development 
of single molecular magnet as qubits^. Therefore, the 
spin engineering in magnetic nanostructures, including 
theoretical modeling, experimental demonstration, and 
device design for applications, is demanding multidisci¬ 
plinary backgrounds of knowledge and technology which 
create great challenges for researchers. 

Magnetic or spin logic appears as an appealing alterna¬ 
tive due to its nonvolatile character, which can boost up 
switching on/off, its possibility to reduce the size of the 
element down to the several-atoms scale on/off spin per 
atom instead of one elementary charge per 10^ atoms in 
semiconductors, and speed increase as a secondary size 
effect^. The modern experiments are promising, they 
still move in the micrometer regime^, thus not fully ex¬ 
ploiting the possible quantum nature of molecular mag¬ 
netism. In other hand, small molecules perform logic op¬ 
erations using as input cation concentrations. The latter 
is fast with respect to the logic operation but slow with 
respect to repeatability. Thus a need for magnetic-logic 


devices on the molecular scale emerges. At the same time 
the experimental evidence of laser-driven ultrafast mag- 
netic(spin) manipulation in the antiferromagnetic mate¬ 
rials motivates the design of a cluster with more than one 
spin center which allows for spin manipulation^i^. 

The problem of the gigahertz magnetization(spin) 
switching speed of todays magnetic logic and a mag¬ 
netic memory devices into the terahertz regime under¬ 
lies the entire field of information processingii^. The 
physics of present-day devices imposes serious limitations 
on this technological transformation, so we must invent 
new paradigms based on the quantum spin dynamics in 
the picosecond regime. This challenge could be met by a 
simulation of the quantum spin switching in a picosecond 
pulsed magnetic field^d^. Currently, no easily accessible 
method is available to generate intense sub-picosecond 
magnetic pulses localized at the nanoscale. However in 
the worbii, show that bimetallic nanorings can act as 
nanoscale sources of intense ultrashort magnetic pulses. 
In this instance rely on the enhanced light absorption as¬ 
sociated with the plasmons of metallic rings to generate 
transient thermoelectric currents that in turn produce 
sub-picosecond pulses with magnetic fields as high as a 
few tenths of a Tesla in the vicinity of the rings. The 
ability to generate strong magnetic fields localized on the 
nanoscale is of interest for elucidating spin and magneti¬ 
zation dynamics at sub-picosecond time and nanometer 
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length scalesi^, and it holds great potential for materi¬ 
als characterization, terahertz radiation generation, and 
magnetic recording. 

Magnetic molecules have been intensively studied so 
far because of their potential technological applications 
and for the possibility they offer to investigate fundamen¬ 
tal properties of matter at the nanoscopic scale^^. An 
intensively studied class of magnetic molecules is consti¬ 
tuted by the so called single molecule magnets or molec¬ 
ular nanomagnets. These systems are characterized by a 
slow relaxation of magnetization at low T and give rise 
to magnetic hysteresis. This is one condition for stor¬ 
ing information in a molecule. Therefore, since a SMM 
represents an isolated magnetic entity, it might be con¬ 
sidered as the smallest practical unit for magnetic memo¬ 
ries. The studying of a possibility of the switching of the 
spin states would incredibly increase the amount of in¬ 
formation storable with respect to current devices. In 
order to have a deeper understanding of fundamental 
features and a major control on technological aspects, 
the time-dependent spin dynamics constitutes the key 
point. The micromagnetic model describes appropriately 
the dynamics of the magnetization(spin) in a magnetic 
nanostructures. In this model the magnetization(spin) 
is assumed to be a spatial- and time-dependent contin¬ 
uous function. The magnetization(spin) dynamics are 
described by the Landau-Lifshitz-Gilbert (LLG) equa¬ 
tion including the energy contributions of the anisotropy, 
the exchange interaction, the magnetostatic interaction, 
and the Zeeman energy!^. The microscopic result of 
the quantum modeling is the cluster spin models, which 
have been proved to be a powerful tool for approaching 
the extreme phenomenology of the ultrafast magnetiza- 
tion(spin) dynamics. Ab initio models for the general 
spin Hamiltonian can be efficient in providing insight on 
the local atomic scale values such as the local magnetic 
moment the local anisotropy D, or a, pairwise ex¬ 
change Jij in the magnetic molecule. The parameters 
of the theory, the exchange integrals and the anisotropy, 
are usually fitted to experiments or calculated from static 
density functional theory (DFT)i^. 

Switching between the LS and HS states can be done 
by varying the temperature, applying pressure or by 
light irradiation. In the first and the second case the 
population of the LS and HS states is fully determined 
by thermodynamic equilibrium and such transitions are 
called temperature- and pressure-induced spin crossover 
transitional^. In the case of light-induced spin crossover, 
after initial excitation to d-states or to metal-to-ligand 
charge transfer states, the system undergoes a complex 
cascade of intermediate transitions and eventually gets 
trapped in a non-equilibrium metastable state. The 
whole process is called light-induced excited spin state 
trapping (LIESST). The details of the LIESST process 
are still under debate^. 

We report here the possibility of a commutation of 
the CoOEP SMM between the two states of a spin- 
crossover system by a pulsed magnetic fieldi^. Choos¬ 


ing a cobalt-based spin-crossover system follows the idea 
that cobalt(H) [CoOEP] materials^i as well as iron(HI) 
ones have metal-ligand bonds shorter than those in 
iron(H) materials, leading to faster relaxation dynamics 
between the HS and LS states. The molecule cobalt(H) 
([Ar]3d'^4s^) may be either in the LS state (tititit S = 
1/2) or in the HS state (tj-tittt S = 3/2), both of which 
are paramagnetic. 


II. THEORETICAL ASPECTS 
A. The spin-dynamics simulations 


The dynamic behavior of a spin is determined by the 
equation of motion, which can be derived from the quan¬ 
tum theory with the general spin Hamiltonian Hgpin that 
calculated the spin structure of a magnetic molecule with 
Hamiltonian 


Hspin — “t” ^an T ^ZEE T 


( 1 ) 


The first term Hex is the Heisenberg-Dirac Hamilto¬ 
nian, which represents the isotropic exchange interac¬ 
tion, Han is the exchange Hamiltonian the term due to 
the axial single-ion anisotropy, and Hzee is the interac¬ 
tion between the spin system and the external magnetic 
field. In order to give a theoretical description of a mag¬ 
netic molecule we exploit the irreducible tensor operator 
technique's. For the study of the CoOEP magnetic prop¬ 
erties, it is a good approximation to treat the electrons 
of one Co-atom as a single atomic spin which is interact¬ 
ing with its surroundings. Thus, one can be content with 
Hgpin without the first term Hex- We will consider an 
approximate treatment, restricting ourselves to the sim¬ 
plest, but very useful, approximation, known as mean 
field theory. This ” mean field” HamiltonianSi, which de¬ 
scribes the interaction of the spin S with the external 
magnetic field, given by its flux He//, can be expressed 
as: 


^spin — 0^^e//S. 

The generally shaped formula is 


spin 

elf = -- 


9 (S) 


( 2 ) 

( 3 ) 


Here the effective magnetic field He// is an external 
magnetic field H^-, the anisotropy fields H^^, the ex¬ 
change interaction Hea; and external magnetic pulse field 
Hp„ise(f)- We use the ’’mean field”— approximation for 
He// ^ H/}®/’" with replacement S ^ = 7(8). We 

call ’’mean” the mean field. (Historically it was often 
called the ’’molecular field”). It includes the average ef¬ 
fect of the neighbors but neglects correlations between 
the spin and its neighbors. Since we have the one mag¬ 
netic atom, this approach us quite suitable. Use (2) we 
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obtain that^ 


B. Calculations of spin-Hamiltonian parameters 


(4) 

The effective magnetic field is given by the free 

magnetic energy variational with magnetization: 


H 


mean 

eff 




6F 

SMs^ 


where F is the free energy of the magnetic nanosystem 


F = -NkeTlnZiH^) 
with the partition function 


Z{F[z)= ^ exp[-e^(Ms)/fcT] ^exp[-geMsi/ 2 /fcr] 

(5) 

. Here we have the energy levels efi{Ms) of the spin- 
Hamiltonian Flspin, which stay we diagonalizied Flgpin- 

A spin structure is defined only proceeding from the 
spin model of a molecule. Here we use to calculate a spin 
structure by the ITO method within the generalized spin 
Hamiltonian Hspin- 

Once we have the energy levels of the spin- 
Hamiltonian, we can evaluate a different thermody¬ 
namic properties of the system as the magnetization, 
the magnetic susceptibility, and the magnetic specific 
heat. Since in further researches the anisotropic part 
of a molecule will be only scalar, the magnetic properties 
of the anisotropic system do not depend on the direction 
of the magnetic field. Thus we can consider the external 
magnetic field iJz directed along arbitrary axis z of the 
molecule coordinate frame that is chosen as a spin quan¬ 
tization axis. In this case the energies of the system will 
be €p{Ms)+geMsHz, where €p{Ms) are the eigenvalues of 
the spin-Hamiltonian containing the magnetic exchange 
and the double exchange contributions (index /i runs over 
the energy levels with given total spin protection Ms). 
Using this expression one can evaluate the magnetization 
Ms by the standart thermodinamic definitions: 




an 


an 


( 6 ) 


The effective field can be derived from the free 

energy functional 


H 


mean 

eff 


S{F{Ms,H,) + F{t)) 

SMs 


dF{Ms,H,) 

aMs 


+H: 


pulse 


(7) 

We have derived a general form of the time-dependend 
spin equation for a system of the spins precessing in an 
effective magnetic field with specifying the interactions 
in the magnetic molecule. 


The present research is devoted to study a 
time-dependent behavior of the Co-octaethylporphyrin 
molecule after applying an external magnetic field in a 
picoseconds time span. In order to optimize the molecule 
geometry and determinate the spin-Hamiltonian param¬ 
eters - zero-field splitting (ZFS, D-tensor), g-tensor, the 
exchange parameters , and the single-ion anisotropy 
parameters, we performed Density Functional Theory 
(DFT) calculations with taking into account the spin- 
orbit coupling effects using ADF package^^i^. The calcu¬ 
lations were performed with the Perdew-Burke-Ernzerhof 
(PBE) exchange-correlation functional and the ZORA 
scalar relativistic Hamiltonian, using an all-electron va¬ 
lence triple zeta with a one polarization function (TZP) 
basis set and a ’Good Becke’ grid. The zero-field split¬ 
ting (ZFS) of the Co OEP ground state was found to 
be equal 1.32meV. ZFS is the breaking of degeneracy 
of the ground state that is not described by a standard 
non-relativistic Hamiltonian. As calculated by ADF, the 
ZFS is that exhibited by molecules whose ground state 
characterized by the spin S>l/2 and absence of a spa¬ 
tial degeneracy. This type of ZFS has two contributions, 
second-order spin-orbit coupling and spin-spin coupling. 
In the present implementation only the spin-orbit cou¬ 
pling term is included. In our study the ZFS was cal¬ 
culated in combination with GGA PBE functional. In 
order to calculate ZFS the relativistic scalar Zero Order 
Regular Approximation (ZORA) option^^ was included. 
The ZFS Hamiltonian is given by 

FIzfs = S • D • S (8) 

where S is the vector operator for the effective spin of 
Co atom, and D is the second-order anisotropy tensor. 
After diagonalization in the magnetic axis frames, the 
Hamiltonian becomes 


Hzfs = DSl + E{Sl - SD (9) 


The parameter D is called the axial ZFS parameter and 
E is the rhombic ZES parameter, which vanishes for high 
symmetry molecules. Parameter D has to be negative for 
the magnetic bistability assosiated with SMMs. When a 
molecule is placed in a magnetic field H, the electron’s 
energies will depend on its magnetic moment value Mj. 
To account for the anisotropy of the Zeeman response to 
an applied magnetic field, an ’’effective” Zeeman Hamil¬ 
tonian using a so-called ”g tensor” is used 
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(10) 


The g-tensor was calculated in a spin-orbit coupled spin 
unrestricted relativistic ZORA approach. The degener¬ 
ate perturbation theory was used with the external mag¬ 
netic field as perturbation. 














4 


III. RESULTS AND DISCUSSION 

The most important aspect to the establishment of real 
devices based on SMMs, e.g. high density rewritable 
memories, is to achieve such conditions under that 
externally driven system could be switched between 
two stable LS and HS spin states and would remain 
in each of these states for a certain period of time. 
The LS/HS bi-stability is connected with small energy 
{AEhs-ls=2-133 meV ) magnitudes involved in the 
switching between these two spin states (Figure la). 
Much research is currently focused on the study of photo¬ 
switching mechanism between LS and HS states, changes 
in the geometric and electronic structure of molecules in 
these states, the magnetic properties of SMMs. How¬ 
ever, we found no studies on the time-dependent spin 
dynamics of spin-crossover systems. In this paper we fo¬ 
cus on the behavior of the spin of the molecules under 
the influence of external stimuli and finding the condi¬ 
tions under which the system would be able to stay in 
one of the bistable states arbitrarily long time. Since 
the photo-induced states of SMMs are characterized by 
a short lifetime, we effect on the system being studied 
by a short magnetic pulse of high amplitude. In the ex¬ 
perimental studjd^, it was shown that the application of 
an external magnetic field of 32T to spin-crossover solid 
system Co{H 2 {fsa) 2 en){py )2 takes place an irreversible 
and quasicomplete transition from the HS to the LS state. 
This also means that after applying a perturbation the 
system can remain in HS state for a long time. In this 
paper, we have developed an approach for the study of 
the spin dynamics of spin-crossover systems and apply it 
to one of the intensively studied objects from a number 
of SMMs the CoOEP molecule. The chemical struc¬ 
ture of CoOEP molecule is exhibited in Pig. HKc). FigH 
shows the magnetic susceptibility plotted as the product 
xT versus temperature for the CoOEP molecule. The 
temperature variation of the magnetic susceptibility yT 
showing the abrupt of the thermal spin transition with 
the jump between two spin states occurs at temperature 
of ri/2=40K. 

By irradiating LS molecules with visible light, HS 
molecule can be generated because of the efficient elec¬ 
tronic decay of excited state via Inter System Crossing 
towards matastable HS state. This effect is now referred 
as LIESST (Light Induced Electronic Spin State Trap¬ 
ping) in the literature^^. Once irradiation is stopped, 
the system relaxes to the thermal equilibrium state LS. 
One important LIESST-result is the decoupling of differ¬ 
ent degrees of freedom during the electronic spin state 
trapping and the transfer of excess energy to molecular 
vibration modes. Contrary to the LIESST, here con¬ 
sidered in some detail the quantum dynamics of the HS 
spin state trapping and impacts on it of a short magnetic 
pulse. 

We present here a theoretical approach to the or¬ 
ganic magnetic CoOEP-molecule in terms of numerically 
solvable the time-dependent LLG-model with a general 
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FIG. 1: (Color online) (a)Spin structure of CoOEP molecule 
in a spin crossover system HS-LS states;(b)HS and LS states 
with Zd’ configuration of in a tetrahedral field;(c)The 

molecular structure of CoOEP. 


spin Hamiltonian in the context of the ITO model. In 
practice, these models apply to the quantum spin-state 
switching in the Co-octaethylporphyrin molecular mag¬ 
net with a picosecond pulsed magnetic field and pro¬ 
vide an understanding of a spin-crossover phenomenon 
in these molecule. In the calculations by LLG-method 
was used the next parameters: Tp„;se=0.I7ps (ultrashort 
THz pulse); time of start pulse is 37.4ps; the height of 
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FIG. 3: (Color online) The spin dynamics for the CoOEP 
molecule with magnetic pulse (S=l/2). 


pulse is 36.ST. We solved the LLG equation for two dif¬ 
ferent states S=l/2, t^gCg (LS) and S=3/2, t^gC^g (HS). 
Under the influence of the constant Zeeman magnetic 
field O.IT, after a certain period of time, which is deter¬ 
mined in case of a single molecule by the size of ZFS, the 
time-dependent spin structure is formed. In the time of 
” saturation” of the system we apply the external stimuli 
of a magnetic pulse. The system in the LS state behaves 
classically without the application of the pulse. Under 
the influence of the 36.8T pulse a spin switching on the 
pulse duration of 100/s is happened, after which the sys¬ 
tem returns to the ground state (FiglS]). In the case of 
the HS state of the system the situation is changing dras¬ 
tically. Under the influence of the pulse a spin switch¬ 
ing occurs as well, but the system remains in this state 



time [ps] 


FIG. 5: (Color online) The spin dynamics for CoOEP 
molecule without magnetic pulse (S=3/2). 


for a long time (Figl6]). With increasing magnetic field 
up to 36.8T a spin switch occurs and the Hgff sign is 
changed. In our approach Heff is the (S) derivative 
of the spin-Hamiltonian. In our case, the Hamiltonian 
m has no Hex term, because there is only one transi¬ 
tion metal atom in CoOEP molecule and, accordingly, 
the spin-spin interaction is absent. Thus, the anisotropic 
Han term is only a nonconstant one, which is quadratic 
in (S) @ and, therefore, He// is proportional to He// is 
the (S). Qualitatively, the He// behavior is the same as 
the (S), which is seeing from the Fig (S] and HI 
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molecule in the LS and HS states. We found that the 
spin-crossover CoOEP molecule undergo the spin switch¬ 
ing between LS and HS states after applying the exter¬ 
nal stimuli of a short magnetic pulse of 36.8T. The tem¬ 
perature dependency of magnetic susceptibility shows an 
abrupt thermal spin transition between two spin states 
at 40K. We also found the difference in the spin dynam¬ 
ics for two spin states of CoOEP molecule. Thus, after 
the spin switching under the magnetic pulse the system 
in the LS state reverts to its ground state, whereas the 
system in the HS state after the spin switching remains 
in the excited state for a long time. This behavior gives 
us the possibility to assign such system as an appropri¬ 
ate candidate for the SMM and opens a way to realize a 
controlled spin switching in real devices for spintronics. 


FIG. 6: (Color online) The spin dynamics for the CoOEP 
molecule with magnetic pulse (S=3/2). 

IV. CONCLUDING REMARKS 

In the paper we present the results of the theoreti¬ 
cal study of the spin behavior in the CoOEP molecule 
under influence of high magnetic pulse. Based on the 
LLG equation we described the spin dynamics of the 
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